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ABSTRACT 



CN ' Context. The emission from Herbig-Haro objects and supersonic molecular outflows is understood as cooling radiation behind 

shocks, initiated by a (proto-)stellar wind or jet. Within a given object, one often observes the occurrence of both dissociative (J-type) 
' ^ ' ■ and non-dissociative (C-type) shocks, owing to the collective effects of internally varying shock velocities. 

' Aims. We are aiming at the observational estimation of the relative contribution to the cooling by CO and H2O, as this provides 

' decisive information for the understanding of the oxygen chemistry behind interstellar shock waves. 

• , Methods. The high sensitivity of HIFI, in combination with its high spectral resolution capability, allows us to trace the H2O outflow 

. wings at unprecedented signal-to-noise. From the observation of spectrally resolved H2O and CO lines in the HH52-54 system, both 

(Ph' from space and from ground, we arrive at the spatial and velocity distribution of the molecular outflow gas. Solving the statistical 

equilibrium and non-LTE radiative transfer equations provides us with estimates of the physical parameters of this gas, including the 
cooling rate ratios of the species. The radiative transfer is based on an Accelerated Lambda Iteration code, where we use the fact that 
variable shock strengths, distributed along the front, are naturally implied by a curved surface. 
^ ' Results. Based on observations of CO and H2O spectral lines, we conclude that the emission is confined to the HH54 region. The 

' quantitative analysis of our observations favours a ratio of the CO-to-H20-cooling-rate » 1 . Formally, we derive the ratio A(CO)/ A(o- 

H2O) = 10, which is in good agreement with earlier determination of 7 based on ISO-LWS observations. From the best-fit model 
to the CO emission, we arrive at an H2O abundance close to 1 x 10"^. The line profiles exhibit two components, one of which is 
^ ' triangular and another, which is a superposed, additional feature. This additional feature likely originates from a region smaller than 

' the beam where the ortho-water abundance is smaller than in the quiescent gas. 

, Conclusions. Comparison with recent shock models indicate that a planar shock can not easily explain the observed line strengths 

. and triangular line profiles. We conclude that the geometry can play an important role. Although abundances support a scenario where 

' J-type shocks are present, higher cooling rate ratios than predicted by these type of shocks are derived. 

■ Key words. Stars: formation - Stars: winds, outflows - ISM: Herbig-Haro objects - ISM: jets and outflows - ISM: molecules 
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1. Introduction 11953 ). tracing the gas at highest velocity. iLiseau & Sandelll 



1986h showed that HH-objects and outflows are physically as- 
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Outflows have many times been discovered through observations ^^^.^^^^ j j . ^j^^^ ^^^^ j^^^^ ^j^^ ^^^^ ^^^^^^ 
of Herbig-Haro obiects (see e.g. He rbig,.1950. i951; Haro 19521 ' fj ^ j j & 



H Water is one of the coola nts that is most sens itive to different 

, * Herschel is an ES A space observatory with science instruments pro- ^j^^^,^ chemistry (e.g. iBergin et aUflHl) . Depending on 

vided by European-led Principal Investigator consortia and with impor- ^. • • ^- r ^- ^ u » .1. j i r 

^ , -f. . ^.^ J. , , , „ , the lonisation traction, magnetic field strength and velocity ot 

tant participation trom NASA. , , , , , , , , ,.m > > 

the shock, water abund ances can be elevated to diiierent levels 

Complementary observations were made with: dHoUenbach et alj[l989l) . In J-type shocks (Jump shocks), where 

Odin is a Swedish-led satellite project funded jointly by the the magnetosonic speed is lower than the propagation of the 

Swedish National Space Board (SNSB), the Canadian Space Agency pressure increase, the involved energies generally dissociate H2 

(CS A), the National Technology Agency of Finland (Tekes) and Centre and all molecules with lower binding energies. As such the wa- 

National d'Etude Spatiale (CNES). ter abundance is generally low in J-type shocks, although it may 

The Swedish ESO Submillimetre Telescope (SEST) located at La reform in the post-shock cooling region once pre-shock densi- 

Silla, Chile was funded by the Swedish Research Council (VR) and the ties are sufficiently high. In C-type shocks however, the pre- 

European Southern Observatory It was decommissioned in 2003. shock gas is partially heated due to traversing magnetic waves 

The Atacama Pathfinder Experiment (APEX) is a collaboration from the post -shock gas, aii d molecules can survive the passage 



between the Max-Planck-Institut fur Radioastronomie, the European of the shock ( lDrainel ll98^ where both the magnetic field and 



Southern Observatory and the Onsala Space Observatory. the gas density is compressed. In this type of shock, the activa- 
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tion barrier for neutral-neutral reactions between molecular hy- 
drogen and oxygen is reached, and the water abundance is ex- 
pected to become enhanced. T his can be both due to the e ffect 
of sputtering from dust grains dKaufman & Neufeldl[T996h and 
due t o high temperature chemistry occurring in the shocked re- 
gion jBergin et al.lll998h . After the shock passage, the enhanced 
water abundance persists for a long time in the post-shock gas. 

The atmosphere of the Earth is opaque at the wavelength 
of the lowest rotational transitions of water as well as most 
other higher excited transitions. Thus, it is not until recent 
years, with the use of space based observatories, these transi- 
tions have been observed succe ssfully. Prev ious missions such 
as SWAS (lMelnicketal.1120001) and Odin dNordhet al.ll2003h 
have put constraints on the water abundance and the dynam- 
ics of molecular outflows ( iFrankhn et al.l 120081 : iBierkeli et al.l 
I2009h . None of these missions, however, provided the spatial 
and spectral resolution that is ava ilable with the Herschel Space 
Observatory (iPilbratt et al.ll2010l) . 

HH 54 is a Herbig-Haro object located i n the Chamaeleon I I 
cloud at a distance of approximately 180pc (IWhittet et aLligQTl) . 
The visible objects in HH54 are moving at high Doppler speed 
out of the cloud tow ard the observer, with veloci ties of the or- 
der 10 - 100 kms ' dCaratti o Garatti et al.ll2009l) . The objects 
show a clumpy appearance due to either Rayleigh-Taylor insta- 
bilities in the flow or due to variability in the jet itself. Another 
possibility is patchy overlying dust extinction. The source of 
the jet is not well constrained. This is further discussed in ap- 
pendix |Aj2] where we also present quantitative arguments for 
identifying IRAS 12553-7651 (ISO-Cha II 28) as the HH54 jet- 
driving source and its associated blueshifted CO outflow. No red- 
shifted emission is observed as what would be expected from a 
bipolar jet (see Section O. On the other hand, the extinction is 
relatively low, something that might allow the redshifted gas to 
flow out essentially unhindered into the low density material at 
the rear side of the cloud. It can however not be ruled out that the 
outflow itself is one-sided and asymmetric. Recent simulations 
show that rapidly rotating stars with c omplex magnetic field s can 
be responsible for such type of flows ( iLovelace et al.ll20Tol) . 

HH54 has previously been observed in va rious lines of 
CO, SiO and H2O using Odin and SEST (iBierkeh et al.l 
I2OO9I R. Liseau, unpublished). Several H2O fines and high- 
J CO lines wer e also observ ed with ISO - LWS and published 
in iLiseau et all {l99d) and iNisini et al.l (Il99d) . During the 
Performance Verification Phase of the Heterodyne Instrument 
for the Far-Infrared (HIFI) instrument, aboard the Herschel 
Space Observatory, the CO (10 -9) transition was observe d. As 
part of the WISH keyprogram dvan Dishoeck et aT]l201 ih . also 
the H20(lio-loi) and H20(2i2-loi) transitions were ob- 
served using HIFI and the Photodetector Array Camera and 
Spectrometer (PACS) respectively. We note that the wavelength 
region, covering the H20(2i2 - loi) transition, also can be ob- 
served with HIFI. 

In this paper, we present observations, both from space and 
ground, of HH 52-54 in spectral lines of CO and H2O. We choose 
to observe this region based on the fact that HH54 is free from 
contamination from other objects, spatially confined and re- 
solved in the infrared regime with the instruments used. Using 
results from observations carried out with SEST, APEX, Odin 
and Herschel, we aim at improving our understanding of inter- 
stellar shock waves. The observing modes and the instruments 
that have been used are described in Section |2] while the details 
of the HIFI data reduction can be found in Appendix iBl The ba- 
sic observational results are summarised in Section [3] whereas 
the interpretations are discussed in Section|4] 
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Fig. 1. Upper panel: False colour map of the CO (10-9) in- 
tegrated intensity in the blue line wing (from -1-2.4 kms"' to 
-30 kms '). The positions of HH54 B, HH53 and HH52 are 
indicated with yellow squares. The readout positions for the on- 
the-fly and raster maps are indicated with white dots. 
Middle panel: CO (10-9) map of the integrated intensity ob- 
tained w ith HIFI in the blue li ne wing overlaid on an Ha im- 
age from lCaratti o Garatti et al. 1 (12009) . Contours are from 5.5 to 
33.6 K km s ' in steps of 3.5 K km s ' . 

Lower panel: A zoom of the CO (10-9) integrated intensity 
overlaid on the H20(2i2 - loi) emission obtain ed with PACS. 
The po sitions of HH 54 B, C and K as indicated in lGiannini et al.l 
(l2006h are indicated with yellow squares. 
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Table 1. Molecular line observations carried out with Odin, SEST, APEX and Herschel. 



Telescope 


Molecule 




Frequency 




HPBW 




Date 


'int 








(GHz) 


(K) 


(") 




(YYMMDD) 


(hr) 


SEST 


CO (2-1) 




230.538 


16.6 


23 


0.50 


970811 - 980806 


3 


SEST 


CO (3-2) 




345.796 


33.2 


15 


0.25 


970811 - 980806 


3 


APEX 


CO (4-3) 




461.041 


55.3 


14 


0.60 


070918 


0.5 


Odin 


CO (5-4) 




576.268 


83.0 


118 


0.90 


050502 - 050620 


4 


APEX 


CO (7-6) 




806.652 


155.9 


8 


0.43 


070918 


0.8 


Herschel-HIFI 


CO (10-9) 




1151.985 


304.2 


19 


0.66 


090726 - 100221 


9 


Odin 


HjOd.o- 


loi) 


556.936 


42.4 


126 


0.90 


090609 - 100420 


12 


Herschel-HIFI 


HjOd.o- 


loi) 


556.936 


42.4 


39 


0.76 


100729 


0.05 


Herschel-PACS 


H,0(2p - 


loi) 


1669.905 


79.5 


13 


N/A 


090226 


0.1 



2. Observations 



2.1.2. PACS 



The observations described in this paper were obtained between 
1997 and 2010 with several different facilities. A summary of 
the observations is presented in Table[T]while the Une intensities 
are listed in Table |2] 



2.1. Herschel 



2.1.1. HIFI 



During the Perform ance Verification Phase of HIFI 
dde Graauw et alj 1201 (ih . the CO (10-9) data were obtained 
on 26-27 July 2009 and 21 February 2010. The H2O data were 
obtained 29 July 2010. The 3.5 m Cassegrain telescope has 
a Full Width Half Maximum (FWHM) of 38" at 557 GHz, 
19" at 1152 GHz and 13" at 1670 GHz, respectively. The 
HIFI H20(lio-loi) spectrum presented in this paper was 
obtained in point mode with position switch using band 1 (490 - 
630 GHz). The OFF spectrum is obtained by a single observa- 
tion of a reference point 10' away. The CO (10-9) HIFI maps 
were obtained in two different observing modes using band 5 
(1120 - 1250 GHz). In the dual-beam-switch raster mode, an 
internal chopper mirror is used to obtain an OFF spectrum 3' 
away from the observed position. In the on-the-fly with position 
switch mode, the telescope is scanning the map area back and 
forth. The data were calibrated using the Herschel Interactive 
Processing Environment (HIPE) version 4.2 and 5.0 for the 
CO (1 0-9) and H20(lio- loi) observations, respectively (lOtll 
I2OIOI) . The data reduction of the HIFI maps is described in detail 
in Appendix IB] For the CO (10-9) observation, data from one 
of the spectrometers on-board have been used. The Wide Band 
Spectrometer (WBS) is an acousto-optical spectrometer with 
a 4 GHz frequency coverage. The channel spacing is 500 kHz 
(0.1 kms ' at 1152 GHz and 0.3 kms ' at 557 GHz). For the 
H20(lio-loi) observation, data from the High Resolution 
Spectrometer (HRS) have also been used. The HRS is an 
Auto-Correlator System (ACS) where the resolution can be 
varied from 0.125 - 1.00 MHz. For this observation it was set 
to 0.24 MHz. Observations from the horizontal (H) and the 
vertical (V) polarisations were combined for both observations. 
The spectra were converted to a scale using main beam 
efficiencies, ?7 mb(l 152GHz) = 0.66 and ?7n,b(557GHz) = 0.76 
(IOIbergll2010l) . 



The PACS spectrograph (iPoglitsch et al.ll2010l) is a 5x5 integral 
field unit array consisting of 9'.'4 square spaxels (spatial picture 
elements). The observations were obtained on 26 February 2009 
in line scan mode and cover the 178.8-180.5 fim region, centered 
on the H2O (2i2 - loi) line at 179.5 fim. The blue channel simul- 
taneously covered the 89.4-90.2 yum region, which is featureless 
and not discussed further Two different nod positions, located 
6' from the target in opposite directions, were used to correct 
for the telescopic background. Data were reduced with HIPE 
version 4.0. The fluxes were normalised to the telescopic back- 
ground and subsequently convert ed to an absolute flux based on 
PACS observations of Neptune dLeUouch et al.ll20l6l) . with an 
approximate uncertainty of ~20 % at 180 fim. The spatial res- 
olution at 180 i-im is nearly diffraction-limited (see Table [TJ. In 
well-centered observations of point sources, only about 40 % 
of the light in the system falls within the central spaxel. The 
H20(2i2 - loi) line is spectrally unresolved in the R - 1700 
spectra (Ay ~ 175 km s"'). 

2.2. Odin 

The Odin space observatory carries a 1.1 m Gregorian tele- 
scope and was launche d into space in 2001 (iNordh et al.ll2"00l 
iHialmarson et al.|[2003h . It is located in a polar orbit at 600 km 
altitude. At 557 and 576 GHz, the FWHM is 126" and 118" 
respectively. The spectra were converted to a r,nb scale using a 
main bea m efficiency, 77mh = 0.9, a s measured from Jupiter ob- 
servations ( IHialmarson et alj|2003h . 

The CO (5-4) observations, suffered from frequency drift 
and have, for that reason, been calibrated, using atmospheric 
spectral lines, acquired during the time i ntervals when Odi n 
observed through the Earth's atmosphere dOlberg et alJ 2003h. 
Since the first publication of the CO (5-4) data in lBierkeli et all 
(I2OO9I) . the frequency calibration scheme has improved. Despite 
this, the velocity scale for this particular observation has some 
uncertainties due to in-orbit variations of the local oscillator unit 
frequency. These variations are most likely caused by slight tem- 
perature changes in the spacecraft during each orbit, due to the 
fact that the Earth is located very nearby. On a velocity scale 
these fluctuations correspond to a ~2 km s"' uncertainty. 

A spectrum, showing a ten tative detection of H2 O ( 1 10 - Iqi) 
at 557 GHz was published in lBierkeU et al.1 (l2009l) . Since then, 
however, additional observations toward HH54 have been car- 
ried out in June 2009 and April 2010 for a total on-source time 
of 12 hours. 
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Table 2. Integrated intensities over the line wings and 1 cr 
uncertainties in parentheses. 



Line 


Sotirce 






mb,rms 






(kms-') 


(Kkms-') 


(mK) 


CO (2-1) 


HH 54 


0.9 to -23.3 


34.6 (0.3) 


224 


CO (3-2) 


HH54 


0.9 to -23.3 


54.9(1.9) 


2021 


CO (4-3) 


HH54 


2.4 to -23.3 


104.4 (3.9) 


773 


CO (5-4) 


HH54 


2.4 to -21.8 


9.3 (0.06) 


20 


CO (7-6) 


HH54 


2.4 to -21.8 


128.1 (25.3) 


4741 


CO (10-9) 


HH 52 






284 




HH 53 






211 




HH 54 


2.4 to -21.8 


33.1 (0.3) 


146 


H20(lio-loi)' 


HH 54 


2.4 to -22.8 


9.2 (0.06) 


20 



Notes to the Table: "The velocity bin size when calculating the rms is 
the same as the channel spacing. 'This refers to the spectra obtained 
with HIFI. 



The observing mode for both observations was position 
switching, where the entire telescope is re-orientated to obtain a 
reference spectrum (10' away in June 2009 and 15' away in April 
2010). The spectrometer used is an acousto-optical spectrometer 
(AOS) where the channel spacing is 620 kHz (0.33 kms ' and 
0.32 kms ' at 557 GHz and 576 GHz re spectively). The dat a 
processing and calibration are described bv lOlberg et al.l (l2003h . 

2.3. SEST 

The observations and calibration of the CO (2-1), CO (3-2), 
SiO(2-l), SiO(3-2) and SiO(5-4) data obtained with SE ST 
were already described in detail by iBierkeli et al.l (l2009l) . to 
which we refer the interested reader. 

2.4. APEX 

CO (4-3) and CO (7-6) data were obtained with the 
APEX/FLASH receiver in service mode in July 2006. In 
this projecQ, both CO lines were observed simultaneously in 
a grid map around HH54B centered on oraooo = 12''55™49'5, 
^2000 = -76°56'23". The selected reference position was rela- 
tively nearby at (120", -120"), which resulted in contaminated 
line spectra close to the cloud LSR velocity (see below). The 
telescope pointing was checked by observing the nearby nova 
X Tra (IRAS 15094-6953). The map was spaced by half the 
instrument beam for the CO (4-3) transition which corresponds 
to 7" at 450 GHz. The data reduction was performed in CLASS 
and the spectra were converted to a r^b scale using the main 
beam efficiencies, 0.60 and 0.43 for the CO (4-3 ) and CO (7-6) 
transitions, respectivelly (iGiisten et al.l l2006l) . The channel 
spacing is 61 kHz (0.04 kms ' at 461 GHz and 0.02 kms ' at 
807 GHz). 



3. Results 

The results from our observations are summarised in Table |2l 
where the integrated intensity over the blue line wing is pre- 
sented. 



' ESQ project code: 077.C-4005(A) 




-30 -20 -10 10 



Velocity [km s ^] 

Fig. 2. H20(lio-loi) spectra obtained with HIFI. The 
blue spectral line is the HRS and the black spectral line 
is the WBS. Both observations were centered on HH54 B: 
ffzooo = 12*'55'"50:3, ^2000 = -76°56'23". The dashed Une in- 
dicates the position of the velocity of the cloud. 



3.1. Herschel 

CO (10-9) is only detected in the HH 54 region. No emission is 
detected toward the region of HH 52-53 down to an rms level 
of ~0.2 K (see upper panel of Fig. [1]) and these objects will 
not be discussed further. From the CO (10-9) emission, we esti- 
mate the size of the source to ~27". The H20(lio - loi) line is 
clearly detected toward HH54 (see Fig.|2]i. Simultaneously, the 
NH3(lo-Oo) line at 572 GHz was covered in the upper side band. 
No emission is detected down to an rms level of ~20 mK. For the 
H20(2i2 - loi) line observed with PACS, emission is detected 
in most of the spaxels and the angular extent of the source is no 
larger than ~28". The peak flux in the central spaxel is 9 Jy. In 
the lower panel of Figure [1] the CO (10-9) integrated intensity 
contours are shown overlaid on the H20(2i2 - loi) normalised 
flux in each spaxel. In this figure, each spaxel is presented on a 
square grid. In reality h owever, there is a sma ll misalignment be- 
tween each spaxel (see iPogUtsch et ani2010l Their Fig 10). An 
off'set of ~9" between the CO (10-9) and the H20(2i2 - loi) 
emission peak is also observed, where the peak of the CO emis- 
sion is located in be tween the B and C clumps as identified by 
ISandeUet all (119871) . This offset might be real given a pointing 
accuracy of a few arcseconds for Herschel. Noteworthy is that 
the different clumps in the region show detectable proper mo- 
tion over a time scale of a few yea rs (e.g. lSchwartz et al.. 19841: 
ICaratti o Garatti et al.ll2006l |2009|) . However, flie CO and H2O 
observations with Herschel were obtained over a time span of 
only one year and it is therefore unlikely that proper motion is 
the cause of the observed offset. The H2O ( 1 10 - loi) spectra ob- 
tained with HIFI are presented in Fig. [2] The line is self absorbed 
by the foreground cloud at ^lsr = +2.4 km s ' . 

3.2. Odin 

The Odin H20(lio- loi) observations carried out on 9 June 
2009 and 20 April 2010 confirmed the previously published de- 
tection with an improved signal to noise. It is this dataset that is 
used for the comparison with HIFI data in the present paper 
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Fig. 3. CO (10-9) map obtained with HIFI. The map shows the spectra toward the region close to HH54 B and the map has been 
regridded with map spacing equal to 9'.'3. Offsets are with respect to HH 54 B; aaooo - 12h55™50:3, 62(m = -76°56'23". The velocity 
scale (wlsr) and intensity scale (Tmb) is indicated in the upper right corner of the map, l/lsr = +2.4 km s"' with a dashed line. 



3.3. SEST 

The CO (2-1) and CO (3-2) maps were centered with a slight 
offset with respect to HH 54 B, viz 4'.'8. The spacing in CO (2-1) 
was 25"while the spacing between the observations in CO (3-2) 
were of the order 15", i.e one full beam width. The number of 
positions observed in the CO (2-1) map and the quality of the 
baselines in the CO (3-2) does not allow us to put constraints 
on the source size. SiO was not detected, when averaging all 
spectra together, down to an rms level of 10 mK for SiO (2-1), 
15 mK for SiO (3-2) and 7 mK for SiO (5-4). The SEST data 
are presented in Fig.|4] 



3.4. APEX 

The quality of the FLASH data is not fully satisfactory, as 
they suffer from off-beam contamination near the line centre. 
However, as we here are mainly focussing on the line wings, 
this should affect our conclusions only little, if at all. In these 
CO line maps, a bump feature (See Sec. 14.1b in the line profile 
was detected in some positions (at ulsr = -7 km s"', see below). 
A reliable source size from these maps can however not be deter- 
mined. In both maps, the feature is detected in a velocity range 
spanning over Ay ^ 7 km s"' (see Fig.|4| . 



Table 3. Parameters used in the CO model 



Parameters held constant 

Distance to source 

CO abundance 

LSR velocity 

Velocity profile 

Shell thickness 

Source size 

Microturbulence 

Gas to dust mass ratio 

Emissivity parameter 

Dust frequency dependence 

Free parameters 

Ht density 
Kinetic temperature 



Aource = 180pc 

X{CO) = 8 X 10-^ 
fLSR = 2.4 km s"' 
v(r) = 20 r/R„y_a km s" 
ArR,„ = 0.12 

^source 30 

Vturb = 1.5 km s"' 
M,,JM,,,, = 100 
'^250 ,,m = 25 cm^g^' 

/? = 1 



n(H2) = 1 X 10^ - 1 X lO** cm-' 
rk,„ = 10 - 330 K 
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Fig. 4. Same as Figure^ but for the CO (2-1), CO (3-2), CO (4-3) and CO (7-6) maps obtained with SEST and APEX. The APEX 
maps have been regridded with a map spacing equal to 8" for CO (4-3) and 9" for CO (7-6). Note that the CO (4-3) spectra may 
be contaminated with OFF beam contribution. 



4. Discussion 

4.1. Observed line profiles 

Common to the observed transitions in CO and H2O is that 
only blue-shifted emission is detected. For all transitions, the 
maximum detected velocity in the line wing is of the order of 
-20 kms"'. A bump-like feature at ^lsr ~ -^ kms"' is also 
clearly visible in the observed CO (10-9) and CO (2-1) spectra 
and possibly also in the CO (3-2), CO (5-4) and H2O (1 10 - loi) 
data. In the CO (10-9) map, this feature is more prominent in 
some positions than others, and most likely it is spatially unre- 
solved to Herschel (see Fig [3]). The ulsr - -7 kms ' compo- 
nent is also clearly visible in the CO (4-3) and CO (7-6) spec- 
tra observed with APEX (see Figure|4]i where the beam sizes are 
13" and 8", respectively. Also in these maps, this component 
seems unresolved, hence it likely originates from a region with 
an angular extent that is smaller than the telescope beams. In the 
CO (5-4) data, where the beam size is 118", the bump-like fea- 
ture is barely visible. This is expected if the beam filling factor is 
small. A position-velocity diagram of the CO (10-9) transition 
shows a trend of higher velocities being detected at lower decli- 
nation and closer to the position of HH 54 B (see Figs. [3]and|5]), 
i.e where the H20(2i2 - loi) emission peaks (see Fig.lTJ. This 
is also clearly visible in the right panel of Fig. |5] where the inte- 



grated intensity of the bump is presented together with the inte- 
grated intensity for the underlying triangular profile. The inten- 
sity maximum of the bump seems offset by ~10" to the south 
from the peak of the bulk outflow emission. The uncertainty at- 
tributed to this offset due to the baseline subtraction is of the 
order ~5". Assuming that the apex of the shock is located close 
to the HH54B position, one would also expect the highest ve- 
locities in this region. 



4.2. Interpretation of the emission line data 

To compute the line profiles for the observed emission, we use 
an Accelerated Lambda Iteration (ALl) code (See App.lCt. The 
ALl code we use is a non-LTE, one dimensional code, assum- 
ing spherical geometry where several subshells are used. The 
number of cells and angles used in the ray tracing can also be 
arbitrarily chosen. In this work, a curved geometry is compared 
with a plane parallel slab to interpret the observed spectral lines. 

The formation of the observed spectral lines occurs most 
likely in shocked gas. Using r esults from detailed models of 
C-shocks, Ne ufeld et al.l (l2006l) presented estimates of the ex- 
citation c onditions for HH 54. From the analysis of H2-rotation 
diagrams, iNeufeld et aO (l2006l) determined the presence of two 



6 



p. Bjerkeli et al.: Herschel observations of the Herbig-Haro objects HH 52-54 




I 



1.5 



1 4 



0.5 



I 




-40 



40 30 20 10 



-10 -20 -30 -40 



RA offset ["] 



Fig. 5. Left panel: Position (in DEC) velocity diagram of the HIFI CO (10-9) data showing the observed intensity. The bump-like 
feature can be seen at ulsr = -1 km s"'. The offset in RA is ". wlsr =+2.4 kms ' is indicated with a dashed line. 
Right panel: Blue contours are the integrated intensity when the bump-like feature is subtracted from the spectra. Contours are from 
3.6 to 28.7 K kms ' in steps of 3.1 K kms '. Red contours show the integrated intensity for the bump. Contours are from 2.6 to 
7.5 K km s ' in steps of 0.6 K km s ' . Underlying colors show the total integrated intensity over the observed line. The positions of 
HH 54 B, C and K are indicated with yellow squares. 



different temperature regimes, viz. at 400 K and at 10^ K, re- 
spectively. 

We used their analytical expressions for the column density, 



A^(H2) = 7.9 X 10 



,20 



n(H2) 



105 



0.5 



\1000/ 



-0.555 



cm 



-2 



and for the velocity gradient, 

nO.5 



^ = 1.6x10^ 
az 



«(H2) 



105 



7- 



\1000 



1.30 



km s ' pc ' . 



(1) 



(2) 



to compute the emission in the CO and H2O lines. Here n(H2) 
is the pre-shock den sity and the other qua ntities have their usual 
meamng. Following ' Neufeldet al.l (l2006l) we take the compres- 
sion factor of 1 .5, which was used by them for H2, but we assume 
that holds for CO as well. In addition, they assume that the frac- 
tional beam filling of the two-temperature components is given 
by the ratio of the column density derived from the rotation dia- 
gram and that given by Eq. [1] 

For these two temperatures, and for pre-shock densitite^ of 
10"* cm and 10^ cm"^, we compute from, Eqs. ([1]) and (|2| in 
slab geometry, the emission in CO and H2O lines, including the 
shapes of the lines. Other parameters (such as CO abundance, 
source size, microturbulence etc.) are given in Table [3] For the 
H2O abundance we assume X(H20) = 1.0x10 '^, consistent 
with th e upper Umit of < 2 x 10 ^ determined bv lNeufeld et alj 
(l2006h . 

The results are presented in Figs. [T] and |9] As seen in the 
latter figure, the computed line strength of the ground state line 
of 0-H2O for a temperature of 400 K and pre-shock density of 
10^ cm is not far from what HIFI has observed. On the other 
hand, such parameters are not in agreement with the observed 
CO (10-9) line, the intensity of which is severely over-predicted. 
As expected, the low-J lines CO (2-1) and CO (3-2) are not very 



^ These were the pre-shock densities used bv lNeufeld et alj j2006h to 
compute the H2 abundance. 



to observer 



e = 30" 




Fig. 6. A cut through the shells of the spherical model described 
in the text. Essentially all radiation originating from the rear 
side of the sphere (grey) is blocked out by the central blackbody 
source. 



sensitive to the changes in temperature from 400 to 1000 K. In all 
cases the computed rectangular line shape is different from the 
observed profiles, which have a pronounced triangular shape. 

4.3. Emission from a curved geometry 

The spatial distribution of t he gas does also infl uence the ob- 
served line profiles (see e.g. iHartigan et al.|[T987l) . For that rea- 
son we also investigate a scenario where the emission originates 
from a curved geometry. To implement this we use a simple 
model that mimics an expanding shell with a diameter of 30" lo- 
cated at a distance of 1 80 pc (see Fig|6]l. The interior of the spher- 



7 



p. Bjerkeli et al.: Herschel observations of the Herbig-Haro objects HH 52-54 



SEST CO(2-1) 



SEST CO(3-2) 







I 






















(.. 







































30 -20 



-10 10 
Velocity [km s~^] 

APEX CO(4-3) 



20 



30 



15 



10 



w 5 



30 -20 



-10 10 
Velocity [km s"^] 

Odin CO(5-4) 



20 



30 



1 ■ 




















8 










6 










4- 
















OP— "-^ 











Velocity [km s ] 




Velocity [km s ] 
APEX CO(7-6) 




S 1 



-10 10 
Velocity [km s"'] 

Herschel/HIFI CO(10-9) 



-20 



-10 10 
Velocity [km s"'] 



20 



30 



Fig. 7. The red (1000 K) and green (400 K) lines represent the model described in lNeufeld et alj ( |2006|) . using pre-shock densities 
of lO'* cm"^ (dashed) and 10^ cm"^^ (solid). The observed spectra from SEST, APEX, Odin and HIFI are plotted with black solid 
lines. Note that the modelled spectra, not have been baselines subtracted. 



ical shell is empty, i.e. at the temperature of cosmic background 
radiation field. This cold sphere occupies 88% of the radius of 
the sphere, based on a shell thickness of 5 x 10'^ cm. This value 
is in between the slab thickness estimated bv lLiseau et al ] (ll996h 
and the analytical expression for the slab thickness at 180 K. 
The shell t hickness is also consistent w ith the cooling length es- 
timated bv lKaufman & Neufeldl (Il996l) . Using this method, we 
block out essentially all the radiation originating from the oppo- 
site side of the sphere. 



To compare with the observed line profiles, spectra are com- 
puted viewing the curved surface from the front. This is sup- 
ported by the absence of detectable SiO emission which is of- 
ten observe d in molecular outfl ows with relatively high veloci- 
ties (see e.g. Nisini et alJl2007h . Shock modelling carried out by 
iGusdorfet all (120081) suggests that sputtering is not very efficient 
in the velocity regime below 25 km s ' . Therefore, assuming that 
the outflow in HH 54 is observed from the front, and that it has 
a small inclination angle with respect to the line of sight, the 
maximum velocity of the molecular gas is likely lower than this. 



Inspection of the CO and H2O spectra show a maximum ra- 
dial velocity of ~20 km s~'. This is also co nsistent with the mod- 
elling carried out bv lGiannini et al. who suggested a C+J 
type shock with a maximum velocity of 18 kms '. The velocity 
in the shell increases linearly from to 20 kms ' (see TableO. 
The true velocity profile is most likely more complicated. In a 
bow shock, the velocity component perpendicular to the jet di- 
rection is expected to be smaller than the component parallel to 
the jet direction. This makes the model somewhat simpler than 
reality. The velocity field within the shocked region probably 
also has a more complicated profile. Furthermore, we assume 
that the emission in all lines stem from the same region of size 
30". The bump-like feature discussed in Section 143] indicating 
a deviation from a linear velocity profile, has not been consid- 
ered in the modelling presented here. The parameters used in the 
model are summarised in Table [3] 

We set up a grid where we vary the H2 density from 10^ 
to 10** cm"-' and the kinetic temperature of the gas from 30 
to 330 K. Thus, our model is steady state and in equilibrium. 
We choose this approach in order to keep the number of free 
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Fig. 8. The six spectra obtained with SEST, APEX, Odin and Herschel compared to the best-fit model (see Sec. 14. 3.11 . 



parameters as small as possible. Also, it is worth noting, that 
we have not achieved a better fit to the bump-like feature (see 
Sec. 14. II) when varying the density and temperature profiles over 
the shells. To find the best fit density and kinetic temperature, 
the reduced x' is minimised, where the difference between the 
observed and the modelled intensity is evaluated in each veloc- 
ity bin (see Fig. [TO]). The CO (2-1), CO (3-2), CO (5-4) and 
CO (10-9) lines are included in the;if^ - minimisation whereas 
the CO (4-3) and CO (7-6) lines are exluded due to the con- 
tamination from the off position and the high noise level of the 
CO (7-6) observation. For all the CO spectra, we only take the 
line wings into consideration, due to the fact that emission from 
the surrounding cloud is clearly visible in the CO (2-1) spectra. 



4.3.1 . Density, temperature and water abundance 

From the curved geometry model, the best fit gas density and ki- 
netic gas temperature are n(H2) = 9 x 10"* cm"-' and r^in = 180 K 
(formally 177 K) respectively. This implies a total H2 mass of 
~1 X 10"^ Mq. The value of the reduced is 2.3. The CO spec- 
tra obtained with Odin, SEST and HIFI are plotted in Figure [8] 
together with the modelled spectra. The model fits the observa- 



tions well, and we conclude that the geometry of the region can 
be a crucial parameter determining the shape of the line profiles. 
The CO (10-9) line observed with Herschel, however, shows a 
slightly more complicated profile than predicted (see discussion 
in Section l4n ). Thus, the triangular form is distorted by the pres- 
ence of the wlsr - -7 kms"' feature. A minor change in the 
kinetic temperature (i.e to 170 K) provides a better fit to the un- 
derlying triangular shape of the CO (10-9) line without affecting 
the lower-J CO lines by much. The computed maximum opti- 
cal depths are <7 x 10"^ for all the modelled CO lines. Using 
the gas density and kinetic temperature obtained from the CO 
modelling, the observed 0-H2O ground state transition is best 
fit with an ortho-watei abundance with respect to H2 of X{o- 
H2O) = 1 x 10"^ (see Figure [TTTi. The observed total flux in the 
map obtained with PACS corresponds to an integrated line in- 
tensity of 20 K kms"'. This is in agreement with the integrated 
intensity from the predicted line profile to within a factor of 2. 



4.3.2. H2O Line profile predictions 

The line profiles for the six lowest rotational transitions of or- 
tho-water have been computed. The abundance is set to 1 x 10"^ 
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Table 4. The predicted integrated intensities, maximum intensities, continuum levels and optical depths for the six lowest rota- 
tional transitions of orf/zo-water. The maximum and integrated intensities are for the baselines subtracted spectra. The foreground 
absorption, visible in the observed H2O (1 lo - loi) spectrum, is not considered in this table. 
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Fig. 9. The red (1 000 K) and green (4() K) lines represent the 
model described in lNeufeld et al.l (l2006h . using pre-shock densi- 
ties of 10"* cm (dashed) and 10^ cm"^ (solid). The observed 
spectrum for the H20(lio- loi) transition with HIFI is plot- 
ted with a black solid line. The 0-H2O abundance is set to X(p- 
H2O) = 1 X 10-^ 



and the lines, that are predicted to be strong enough to be read- 
ily detected with HIFI, are displayed in Figure [12] Also in this 
figure, red-shifted emission originating from the opposite side 
of the sphere is present for the model with a curved geometry. 
Significant changes in the excitation temperature in the inner and 
outer shells show up in these spectra as weak absorption features. 
For this reason, the spectra have been computed using more than 
100 shells to avoid any drastic changes in excitation temperature, 
due to optical depth effects between each shell. Simultaneously 
with the CO (10-9) observation, H20(3i2 - 221) at 1153 GHz 
was also observed. The observed noise level of ~0. 1 K is how- 
ever too high to detect the line with a predicted strength of 
only ~10 mK. A summary of these predictions is presented in 
Table |4] In this table, the observed integrated intensity in the 
H20(lio - loi) line is 30 % lower than the predicted value due 
to the absorption from foreground gas. If this gas is at a low tem- 
perature, however, the higher transitions should not be affected 
by much. 



5.1 




140 160 180 200 220 
Temperature [K] 



Fig. 10. Contour plot showing the sum of the chi-squares, com- 
paring the modelled spectra and the observations. The 5%, 10%, 
20% and 30% deviations from the minimum value are indicated 
with solid lines. The best-fit density and temperature are indi- 
cated with a cross. 

4.3.3. Cooling rate ratios 

For the curved model described in Sec. 14.31 we derive the cool- 
ing ratio A(CO)/A(o-H20) = 10 where the 557 GHz line is the 
dominant contributor to the water cooling. This value agrees well 
with the earlier determination of 7 based on ISO-LWS data and 
presented bv iLiseau et al., (,1996i) . 

4.3.4. The bump-like feature 

As already discussed in Sec. 14.11 the bump-like feature seems 
unresolved to Herschel, i.e the size of the source is uncer- 
tain. A source size comparable to the beam size of Herschel, 
at 1152 GHz (19"), would have to have a high temperature 
(>1000 K) and low density (~10"' cm"^) to fit the observations. 
This implies that the flux in the high-J CO lines would be highe r 
than what was observed with ISO-LWS (iGiannini et al.ll2006l) . 
However, a source size as small as ~\" , where the temperature 
and density would have to be ~400 K and ~10^ cm"-^ respec- 
tively, is a plausible scenario. In that case however, the ortho- 
water abundance has to be less than 10"^ to fit with the ISO- 
LWS observations. In addition, the source is likely not smaller 
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Fig. 11. The HIFI H20(lio - loi) spectrum compared with the 
model spectrum for X{o-H20) = 1 x 10"^. Spectra with the o- 
H2O abundance increased by a factor of two (dashed line) and 
decreased by a factor of two (dotted line) are also included in- 
cluded. Note that the foreground gas is visible in absorption in 
the observed spectrum. 
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Fig. 12. Line profile predictions for the H2O (2i2 - loi), 
H2O (221 - 2i2), H2O (3o3 - 2i2) and H2O (3i2 - 3o3) spectra as 
observed with HIFI. 



than 1". Also in this case, the high- J CO lines would be stronger 
than what is actually observed. 

Using a source si ze of 10" (th i s size is estimated from 
H2 maps presented in iNeufeld et alj (l2006h ). the observed CO 
"bullet"-emission is well fit with a temperature, Tkin - 600 K, 
and a density, «(H2) = 2x10"^ cm This yields a column den- 
sity, A^(H2) = 5 X 10^° cm"^. For an ortho-water abundance of 
X(o-U20) - 1 X 10"^, the H20(3o3 -2 12) emission is entirely 
accounted for by the ISO-LWS observations. In this case the 
H20(2i2 - loi) emission would be merely about 25% of what 
was observed with PACS and no significant contribution would 
be observed in the H2O (1 10 - loi) line- Therefore, we conclude 
that the 0-H2O abundance in the "bullet" is likely lower than 




Rot. quant, number [J] 

Fig. 13. Observed CO line fluxes (squares) as a function of the 
rotational quantum number J. The errorbars are the estimated 
calibration uncertainties (10% for SEST, Odin and HIFI and 
30% for APEX and IS O). T he best fit model, described in the 
present paper (see Sec. 14.3.11 . is indicated with a solid line. The 
fit to the high-J CO lines is indicated with black circles (see 
Sec. l4331 l. 

4.3.5. Comparison with other results 

The derived H2 density implies a column density of 
5 X 10^° cm"^ which is approximately one order of magnitude 
highe r than what is derived for the warm gas (iNeufeld et al.l 
l2006h . The best fit temperature, - 180 K, is significantly 
lower than the temperature in the gas responsible for the high- 
J CO emission [i.e CO(14-13) - CO (20-19)] reported in 
iGiannini et alj (|2006|) . Our model predicts these lines to be one 
to two orders of magnitude weaker. Clearly, different tempera- 
ture regimes are present in HH 54 and a one-temperature, one- 
density model cannot explain all the infrared observations. The 
high-J CO emission observed with ISO can be explained using 
the same geometry and source size but having a shell thickness 
4 X 10'"^ cm, i.e one order of magnitude thinner. In that case a 
slightly higher density («(H2) = 1.5 x 10^ cm"^) and tempera- 
ture (T = 500 K) fits the high-J CO lines well. On the other hand, 
this secondary component also makes a significant contribution 
to the CO (10-9) line and this emission may origin ate from dif- 
ferent components. Recently, iTakami et al.l (1201 Ol) report mor- 
phological differences between Spitzer observations in the 3.6, 
4.5, 5.8 and 8.0 /im bands. The emission is observed to be less 
patchy in the long wavelength bands (i.e. 5.8 and 8.0 yum) and 
they interpret this as thermal H2 emission being more enhanced 
in regions of lower density and temperature. Given the fact that 
hot gas obviously is present in this region, the secondary com- 
ponent may in reality be in smaller regions of high temperature 
similar to those observed in H2. In Figure[T3] the CO line flux is 
plotted as a function of the rotational quantum number, J. 

The predicted integrated intensiti es can be co mpared with 
the modelling presented in Gianni ni et al.l (l2006h . These au- 
thors present a multi-species analysis where they conclude that 
the observed H2, CO and H2O lines can only be explained 
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by a J-shock with magnetic precursor. The C+J shock model, 
presented in their paper, explains the observed H20(2i2 - loi) 
and H2O (3o3 - 2i2) line fluxes of 7 x lO"'-' erg cm"- s"' and 
2 X 10 '"^ erg cm"^ s ' well, and they predict the line flux for the 
H20(lio - loi) line to be 2.0 x 10"'^ erg cm"^ s"'. Converted 
to a K kms ' scale these values correspond to 6.8, 1.4 and 
53 K kms"' respectively. Taking the beam size into account, 
the predicted integrated intensity for the H2O (1 10 - loi) line in 
this paper is at least a factor of four lower This could be due 
to the fact that this line is very sensitive to the type of shock 
presen t. This is also discussed in the paper by iGiannini et al.l 
(I2OO6I) . where they note that a J type shock would change the 
flux by more than a factor of two downwards. 

Recently, iFlower & Pineau des ForetsI (1201 Ol) presented the- 
oretical predictions of CO and H2O line intensities, based on 
detailed C- and J-type shock model calculations. Expectedly, 
the assumed magnetic fields are different for their models of C- 
and J-shocks, i.e. for b - I and b = 0.1, respectively, where 
b is defined through B - b y/nn jiG. For their high density C- 
shocks (see below), this means that pre-shock fields of order 
450 should be present. Based on OH Zeeman observations, 
iTroland & Crutche j ( l2008h detected 9 dark clouds in a sample of 
34. Corresponding line-of-sight magnetic field strengths were in 
the range Z^io s - 10-25 ^G, with typical v alues around 15 jiG. As 
discussed bv lTroland & Crutcheij (l2008h . for randomly oriented 
fields, Bios - 0.5x|B|. Hence, magnetic fields in dark clouds 
do not likely exceed levels of 20 - 50 ^G. This refers to the ob- 
served scales of 3'. However, these authors also showed that in 
a given cloud, Bios did not change appreciably from one posi- 
tion (active molecular outflow) to another (quiescent surround- 
ing cloud). It seems, therefore, that fields as strong as 450juG 
may not be that common. On the other hand, order of magnitude 
lower field strengths (45 /^G for b - 0.1) are more consistent 
with the observational evidence and may promote the occurrence 
of J-shocks. 

For the radiative transfer lFlower & Pineau des ForetsI (I2OIOI) 
used an LVG approximation in slab geometry. In their Fig. 8, 
line profiles for CO (5-4) and H2O (1 10 - loi) are shown for C- 
shocks with two densities and four shock velocities. These lines 
are close in frequency and observations of HH 54 with Odin are 
made with essentially the same telescope beam, rendering reso- 
lution issues to be of only minor importance, i.e. any scaling of 
the intensities due to the source size should affect both observed 
lines in the same way. It should be feasible, therefore, to directly 
compare the line profiles of our Odin observations with those 
of the models bylH ower & Pineau des ForetsI ( |2010| ). Based on 
the observed maximum radial velocities, we consider only mod- 
els with Ws > 20 kms"'. Their models with shock velocities of 
20 kms"' could correspond to a head-on view, whereas their 
models for = 30 and 40 km s"' could correspond to inclina- 
tions of the flow with respect to the line of sight of 48° and 60°, 
respectivel}0. 

For the C-shock models, several of the line profiles display 
shapes that are qualitatively similar to those observed with Odin 
and HlFl. These line shapes are a consequence of the computed 
flow variables, not a geometry effect. However, in particular for 
the range of Aw ~ to about 10 kms"', the Odin/HIFI lines 
exhibit CO-to-H20 intensity ratios very much in excess of unity, 
i.e. 7'co(Ai')/rH,o(Ay) » 1. Au is the velocity relative to the 
rest frame of the flow, i.e. wlsr- 

' For the emission knots of HH 54 lCaratti o Garatti et al.l l i2006l) de- 
termined an average inclination of 27°, which would imply that ~ 




Velocity [km s ^] 



Fig. 14. Comparison between the CO (5-4) line (black) and the 
H20(lio - loi) line obtained with Odin (red). The H2O spec- 
trum has been multiplied with 6 for clarity. Also the HlFl 
H2O (1 10 - lot) spectrum has been plotted in this figure for com- 
parison (blue). The intensity of the emission has in this case been 
multiplied by 6 and corrected for beam filling. 




Fig. 15. Line ratios between CO (5-4) and H2O (1 10 - loi)- The 
ratio as measured from Odin is indicated with a black fine. The 
predictions presented bv lFlower & Pineau des ForetsI (l2010l ) for 
i^shock =20 kms"' are plotted in blue for «(H) = 2x 10^ cm"^ 
and in red for n(H) = 2x10^ cm"^. 



These observed ratios are very much larger than the theoreti- 
cal values (Fig.[T4land[T5Tl. Albeit the ratios are larger than unity 
for the low-density (nu = 2x10'^ cm"^) C-shocks, these fall still 
far below the observed ones. On the other hand, the high den- 
sity («H = 2 X 10^ cm"-') cases could directly be dismissed, as 
these tend to show inverted ratios, i.e. TcolTn^o «: 1, contrary 
to what is observed (Fig.fTTt. 

For the J-shock models, iFlower & Pineau des ForetsI (1201 Ol) 
list the predicted integrated intensities. Also in this case, the pre- 
dicted CO-to-H20 ratios (for a shock velocity of 20 kms"') are 
much lower than the observed ratio, viz 0.06 and 0.005 for the 
pre-shock densities riu - 2 x 10'*cm"-'and «h = 2 x 10^ cm"-' 
respectively. 
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4.4. Implications for future work 

The modelling of planar C-shocks show that low densities are re- 
quired for cooling rate ratios, A(CO)/A(o-H20) > 1, in contrast 
to our own findings, where densities at least as high as 10'' cm"'' 
seemed implied by the observations (Sect. l43?T] ). The cause for 
this mismatch is not clear to us, but one of the reasons could be 
the difference between modelled planar and curved geometry. 
Such modelling should therefore be attempted. The relatively 
high cooling rate ratio, A(CO)/A(o-H20) = 10, is also not 
easily reconcilable with the presence of a C-type shock where a 
ratio, A(CO)/A(o-H20) <s; 1 is favoured. On the other hand, J- 
shock models predict ratios even smaller than for C-shock mod- 
els. However, the relatively low ortho-water abundance deter- 
mined from the modeling indicate that J-type shocks may con- 
tribute. 

As already discussed in Sect. 14.3.51 different temperature 
regimes are present in the HH 54 region. It would be adequate 
therefore, to repeat the ISO-LWS observations, using the higher 
sensitivity and spatial resolution provided by PACS. The mid- J 
CO lines fall in the wavelength range covered by the Spectral 
and Photometric Imaging Receiver (SPIRE) and these should be 
observed. 



5. Conclusions 

Based on spectral mapping with Herschel of the region contain- 
ing the Herbig-Haro objects HH52 to HH54 we conclude the 
following: 

• The CO (10 - 9) 1152 GHz Hne was clearly detected only 
toward the position of HH 54 with a FWHM i 30", compa- 
rable to the extent of the HH-emission knots in the visible 
and infrared. 

• The H20(2i2 - loi) 1669GHz line was clearly detected to- 
ward HH 54 and with a similar extent as the CO (10 - 9) line. 
An offset of 9" (2.4 x 10'^ cm) between the two species is 
observed, but the reality of this can at present not be firmly 
assessed. The H20(lio - loi)557 GHz line was also clearly 
detected toward HH 54. 

• The CO (10 - 9) spectra show only blueshifted emission, 
with maximum relative velocities in excess of -20kms"'. 
The line profiles exhibit typically a triangular shape, which 
in certain positions is however contaminated by a bump-like 
feature at i/lsr = -7 km s"' . This feature is constant in veloc- 
ity and width. It is limited in spatial extent and may be iden- 
tified by what is commonly called a "bullet". Comparison of 
the observed spectra with analytical bow shock line profiles 
limits the viewing angle to greater than zero but i 30°. 

• The bump is clearly seen in position-velocity cuts, reveal- 
ing two peaks, in addition to a smooth velocity gradient of 
about lO^kms"' pc"'. The low- velocity peak appears close 
to the ambient cloud velocity, whereas the second peak cor- 
responds to the bump. For this feature we determine from 
Gaussian fit measurement of the CO (10 - 9) data a relative 
peak-r^b = 1 .5 K, ulsr = -7 km s"' and FWHM = 5 km s"' . 

• These line features are also observed in the CO (5 - 4) spec- 
trum observed with Odin and in spectra obtained from the 
ground, i.e. in maps of the (2-1), (3-2), (4-3) and (7-6) CO 
transitions. In addition to the triangular line shape, these data 
do therefore also confirm the reality of the spectral bump. 

• We initially use physical parameters for shock models of 
HH54 found in the literature to compute the CO spectra. 
These models presented two-component fits, given by ana- 



lytical expressions for temperature and density in the respec- 
tive slabs. These models were only moderately successful in 
reproducing the observational results, in particular what re- 
gards the line shapes. Considerable improvement was found, 
however, in computed spectra using a curved geometry in- 
stead. 

• Using the best fit model parameters (in a ;t'^-sense) for the 
CO data we computed spectra for H2O. This model fits the 
H2O557GHZ line observed with HIFI for an ortho- water 
abundance with respect to molecular hydrogen, X{H20) = 
1 X 10 ^ 

• A cooling rate ratio, A(CO)/A(o-H20) » 1, is not eas- 
ily reconcilable with recent shock modelling. On the other 
hand, a relatively low water abundance (~10"^), supports a 
scenario where J-shocks contribute significantly to the ob- 
served emission. This is also consistent with magnetic field 
strengths measurements toward dark clouds, where Z?-values 
lower than what is needed for C-type shocks, typically are 
observed. 

• Comparison of our Odin data with line profiles from recent 
detailed C-shock model computations in the literature sug- 
gests that some refinements of these models may be required. 
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A.1. Extinction by dust and dust parameters 

At the position of HH 54, the molecular cloud is relatively ten- 
uous as the dust extinction through the cloud amounts merely 
to Av ~ 2-3 mag {Gregorio-Hetem et al. 1989; Cambresy 199^). 
With the relations of Eqs (2) and (8) of iHavakawa et al.1 (l200ll) . 
this results in a total Ha-column de nsity of A^(H2) =2x 10^' cm"^. 

From [Fe 11] line observations. lGred^ (1 19941) determined the 
visual extinction toward the HH object as Ay -1-3 mag, which 
could suggest that the blueshifted HH 54 is located at the rear 
side of the cloud, a circumstance that potentially explains the ap- 
parent absence of redshifted emissioifl However, an T^y-value of 
> 5, which is significantly higher than that used bv lGredel(ll994l) 
and by others, needs to be invoked to d escribe the Cha II dust 
(see the discussion bv lAlcala et al.l2008l and references therein). 
This reflects a predominantly 'big-grain' contribution to the ex- 
tinction (almost 2 magnitudes larger Ay per unit /iB-v) and we 
consequently chose a shallower dust opacity dependence on the 
wavelength in the FIR and submm, identified by the parameter 
/3 in Table 3, i.e. /3 - I as opposed to /3 - 2 for the gene r al ISM 
For the normalisation wavelength used by iHildebrandl ( Il983h . 
i.e. Ao = 250 urn, we adopt the grain opacity kq - 25cm^g"' 
(lOssenkopf & Henning 1994). The corresponding mass absorp- 
tion coefficient k at long wavelengths is then readily found from 
^,1 = Ko{Ao/Af. 

Combining our own results with those found in the literature 
yields an SED of a putative core at 10 K which provides a strict 
upper limit to the luminosity. Any associated central point source 
would have a limiting luminosity of Lboi < 0.4 Lq. Assuming 
optically thin emission at long wavelengths and using standard 
techniques, the mass of the 10 K core/envelope is estimated at a 
mere 2 Mg, excluding the possibility of a very young age of the 
hypothetical central source. 



A.2. Tlie exciting source 

Several IRAS point sources are found in the surroundings of 
HH 54 and it is plausible to assume that one of these would be the 
exciting source of the HH object. Most often, the exciting source 
is a young (f Myr) stellar object, situated at the geometrical cen- 
tre of the bipolar outflow and this seems also to be the case for 
the CO outfl ow associabl e with IRAS 125 15 - 764 10 near HH 52 
and HH 53 (lKneelll992h . For HH54, the situation is diff'erent 



however: |^ed (ll992|) identified IRAS 12522 - 7640, which is 



positionally coincident with HH 54, as the driving source of a 
collimated blueshifted outflow. It seems however very unlikely 
that this IRAS source is (proto-)stellar in nature. It was de- 
tected essentially only, and barely, in the 60 //m-band, and it was 
shown that thi s IRAS flux is enti rely attributable to [O i] 63 fim 
line emission (iLiseau et al.lll996l) . More recent observations by 
Spitzer are in support of this conclusion , as these were una ble to 
reveal any point source at this position (I Alcala et al.ll2008l) . 

On the basis of proper motion measurements and morphol- 
ogy considerations, Caratti o Garatti et al. (2009) proposed the 
Class I object IRAS 12500 - 7658 to be the exciting source of 
HH54. The IRAS source is located 20' south of HH54, i.e. at 
the projected distance of 1 pc. It was not detected at 1200 //m 



From H2-line emissi on at various positions of the HH ob- 
ject, I Giannini et alj i2006l) derived the essentially constant value of 
Av = 1 mag. 

5 lAlcala et al.l ( l2008h identify this source as a likely background K- 
giant. 
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by lYoung et alJ (I2005L see below). What regards HH52 and 
HH53, no satisfactory candidate was found by these authors. 
The po ssibility of differen t ex citing sources was also consid- 
ered bv lNisini etaP d 19961) and lAlcala et al.l (l2008h . associating 
HH52 and HH 53 wit h IRAS 12496 - 7650 (DKCha), whereas 
lib arra & Ladal (|2009|) identify this object also with the exciting 
source of HH 54. 

In APEX-LABOCA maps at 870jum, no emission was de- 
tected toward HH 54 nor from the Cha II cloud itself (S svo/jm < 
20mJy beam Ivan KempenI l2008l) . However, jet-like exten- 
sions are seen to emanate from both DKCha and IRAS 12553 - 
7651 (ISO-ChaII28), which are situated some 10' to 20' 
south/southeast of HH 54. Narrow, ext ended emissi o n from these 
IRAS sources was also reported by lYoung et al.l (l2005h . who 
used SEST-SIMBA to map the region at 1200 yum. The align- 
ment with the HH objects is poor, however Based on the ob- 
served spectral indices, this emission is dominated by optically 
thin radiation from dust. The dust features appear as tori or rings, 
seen edge-on, but given their sizes (3x10"* AU), they do not rep- 
resent what commonly is called a (protostellar) "disc". Perhaps, 
these dust features are left-overs from an earlier history, during 
which the IRAS sources were formed. A "conventional disc" of 
size ~ 10^ AU could very well hide inside and any outflow would 
be orthogonal to these disc features. 

Whereas the dust extension from DKCha is not, whether 
direct or perpendicular, pointing anywhere near HH 52-54, the 
vector orthogonal to the linear feature of ISO-ChaII28, at po- 
sition angle 51°, is within 5° from the current direction toward 
HH 54. Taking the uncertainty of this estimate into account, this 
coincidence appears very compelling. If this source drives/has 
driven a jet, that would very well be aligned with HH 54. The 
projected distance is 16' (0.8 pc), i.e. jet travel times would be 
of the order of 10"* (i//100kms"')"' yr. With a total luminosity 
of ~ 10 Lq arid a m ass accretion rate of Mace = 6 x 10"^ Mq yr"' 
(lAlcala et al.l2008h . ISO-Cha II 28 would make an excellent can- 
didat^ for the exciting source of HH54, the (distance cor- 
rected) stellar mass loss rate for whic h has been determined as 
Mioss < 4 X 10-^ Mq yr-' (lKneelll992h . assuming a flow velocity 
of > 1 00 km s" ' . A ratio of loss-to-accretion rate of < 2/3 would 
thus be indicated, a value which is in reasonable agreement 
with theoretical ex pectat ion for cen trifugally driven MHD winds 
dKonigl & PudritdboOOt fShu et aUlIOOOl: IContopoulos & Sautvl 
120011) . To summarise, as also concluded by others before, the is- 
sue of the exciting source(s) of the HH objects 52, 53, 54 remains 
essentially unsettled. 

Appendix B: HIFI data reduction 

The CO (10-9) data published in this paper are public to the 
scientific community and can be downloaded from the Herschel 
Science Archiv^ (HSA). The two CO (10-9) observations were 
carried out in dual-beam-switch raster mode (observation id: 
1342180798) and on-the-fly with position switch mode (obser- 
vation id: 1342190901). The H20(lio - loi) observation was 
carried out in point mode with position switch. The data reduc- 

* From the data for CO outflows compiled by IWu et"an i2004h . and 
complementing information in the hterature, it is found that Z,(,oi J^ioL' 
where x = 1.0 for Mio^s < lO^^^Moyr"' (for v„-,„i =100 km s"') and 
which steepens to x = 1.2 beyond that. This holds over eight orders of 
magnitude, albeit with considerable spread due to the inhomogeneity 
of the data. According to this relation, Lbd = IOLq if Mioss = 4 x 
10-^Moyr-'. 

' http://herschel.esac.esa.int/Science_Archive.shtml 



-20 20 




10 -10 

RA offset ["] 



Fig. B.l. CO (10-9) spectra toward HH 54. The black spectra are 
from the HIFI observation 1342180798, while the blue spectra 
are from the observation 1342190901. The horisontal and verti- 
cal polarisations have been averaged together (see text). 

tion and production of the HIFI maps followed the following 
steps: 

1 . The raw data were imported to HIPE from the HSA. 

2. All data were re-calibrated using v4.2 of the HIFI pipeline 
for CO (10-9) and v5.0 for H2O (1 lo - loi). There is a sHght 
offset between the pointings for the horisontal and verti- 
cal polarisations. For that reason, separate readout positions 
were calculated for the two polarisations. 

3. The level 2 data were exported to Class using the HiClass 
tool in HIPE. 

4. A baseline was fitted to and subtracted from each individual 
spectrum. This step and the subsequent data reduction was 
carried out in xfl 

5. From Figure [Tj it is clear that the readout positions for 
the two CO (10-9) maps are not on a regular grid in RA 
and Dec. Instead of averaging nearby spectra together, a 
Gaussian weighting procedure was used. The FWHM of the 
two dimensional Gaussian used in the weighting was taken 
to be 9.3". This size is larger than the average spacing be- 
tween individual spectra but on the other hand small enough 
not to change the spatial resolution significantly. 

6. In order to verify the quality of the CO (10-9) observations, 
spectra taken toward the same positions in the two differ- 
ent observing modes were compared, showing no signifi- 
cant differences, (see Figure IB. lb . Also the horisontal and 
vertical polarisation data were reduced individually for both 
CO (10-9) and H20(lio-loi) showing excellent agree- 
ment. 

Appendix C: Radiative transfer analysis 

An Accelerated Lambd a Iteration (ALI) method (e.g. 
iRvbicki & Hummer! 1199 lb is used to compute the line pro- 
files for the observed transitions, varying only the density, 
kinetic temperature and the geometry of the shocked region. 
ALI is a method, where the coupled problem of radiative 
transfer and statistical equilibrium is solved exactly. 

Data reduction software developed by P. 
Bergman at the Onsala Space Observatory, Sweden; 
pittp://www.chalmers.se/rss/oso-en/observations/| 
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C.1. Accelerated Lambda Iteration 

The ALI code has in recent years been used in several pub lica- 
tions (see e.g lJusttanont et al] 120051: fWirstrom et alJ 20Toh and 



was be nchmarked with other codes in a paper bv lMaercker et af] 
( I2OO8I) . In that paper, the ALI technique is described in more de- 
tail. In the present work, typically 30 shells and 16 angles are 
used. The statistical equilibrium equations and the equation of 
radiative transfer are solved iteratively until the relative level 
populations between each iteration changes by less than 10"^. 
In the modelling we set the energy limit to 2000 K. 27 rotational 
energy levels, 26 radiative transitions and 351 collisional tran- 
sitions are included in the calculations for CO. For 0-H2O, the 
corresponding numbers are 45, 164 and 990. 



C.1.1. Collision rates 

In recent years, substantial efforts have been made to deter- 
mine cross-sections for H2O excitation due to collisions with H2. 
Quantal calculations governing collisions between H2O and H2 
are to this date not complete and for 0-H 2O, only collision rate s 
for interaction with P-H2 are available jDubernet et alj l2b09h . 
For this reason we make the assumption that all /7-H2 are in 
the lowest energy state when solving the statistical equilibrium 
equations. It is also assumed that the /?-H2 stays in the ground 
state after the colUsio n. The ortho-to-para ratio has been esti- 
mated for warm H2 bv lNeufeld et al.l (|2006|) . These authors de- 
rive ratios significantly lower than 3, viz ~0.4 - 2 in the HH54 
region. The computed spectra from the ALI modelling have also 
been compared wi th those, using the collision rates presented by 
iFaure eta l] (l2007h . and with similar results. In the case of CO, 
we use the collision rates presented bv lYang et alJ (l2010l) . 



16 



